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(54) Generating and displaying spatially offset sub-frames 



(57) A metliod of displaying an image (1 4) witli a dis- 
play device (26) Includes receiving image data (16) for 
the image. A first sub-frame (30) and a second sub- 
frame (30) corresponding to the image data are gener- 
ated. The first and the second sub-frames are generated 



based on minimization of an error between the image 
data and a simulated image (412, 512, 610, or 706). The 
method includes alternating between displaying the first 
sub-frame in a first position and displaying the second 
sub-frame in a second position spatially offset from the 
first position. 
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Description 

Cross-Reference to Related Applications 

5 [0001] This application is related to U.S. Patent Application Serial No. 10/213,555, filed on August 7, 2002, entitled 
IMAGE DISPLAY SYSTEM AND METHOD; U.S. Patent Application Serial No. 10/242,195, filed on September 11, 
2002, entitled IMAGE DISPLAY SYSTEM AND METHOD; and U.S. Patent Application Serial No. 10/242,645. filed on 
September 11 , 2002, entitled IMAGE DISPLAY SYSTEM AND METHOD. Each of the above U.S. Patent Applications 
is assigned to the assignee of the present invention, and is hereby incorporated by reference herein. 

10 

The Field of the invention 

[0002] The present invention generally relates to display systems, and more particularly to generating and displaying 
spatially offset sub-frames. 

15 

Baclcground of the Invention 

[0003] A conventional system or device for displaying an image, such as a display, projector, or other imaging system, 
produces a displayed image by addressing an array of individual picture elements or pixels an^anged in horizontal rows 
20 and vertical columns. A resolution of the displayed image is defined as the number of horizontal rows and vertical 
columns of individual pixels fomning the displayed Image. The resolution of the displayed Image is affected by a reso- 
lution of the display device Itself as well as a resolution of the Image data processed by the display device, and used 
to produce the displayed Image. 

[0004] Typically, to Increase a resolution of the displayed image, the resolution of the display device as well as the 
25 resolution of the Image data used to produce the displayed image must be increased. Increasing a resolution of the 
display device, however, increases a cost and complexity of the display device. In addition, higher resolution image 
data may not be available and/or may be difficult to generate. 

Summary of the invention 

30 

[0005] Onefonn of the present Invention provides a method of displaying an Image with a display device. The method 
includes receiving image data for the Image. A first sub-frame and a second sub-frame coresponding to the image 
data are generated. The first and the second sub-frames are generated based on minimization of an error between 
the image data and a simulated image. The method includes alternating between displaying the first sub-frame in a 
35 first position and displaying the second sub-frame in a second position spatially offset from the first position. 

Brief Description of the Drawings 

[0006] 

40 

Figure 1 is a blocic diagram illustrating an Image display system according to one embodiment of the present 

invention. 

Figures 2A-2C are schematic diagrams Illustrating the display of two sub-frames according to one embodiment of 
the present invention. 

4s Figures 3A-3E are schematic diagrams Illustrating the display of four sub-frames according to one embodiment of 

the present invention. 

Figures 4A-4E are schematic diagrams illustrating the display of a pixel with an image display system according 
to one embodiment of the present invention. 

Figure 5 Is a diagram Illustrating the generation of tow resolution sub-frames from an original high resolution image 
so using a nearest neighbor algorithm according to one embodiment of the present Invention. 

Figure 6 is a diagram illustrating the generation of low resolution sub-frames from an original high resolution Image 
using a bilinear algorithm according to one embodiment of the present Invention. 

Figure 7 is a block diagram illustrating a system for generating a simulated high resolution Image according to one 
embodiment of the present Invention. 
55 Figure 8 is a block diagram illustrating a system for generating a simulated high resolution Image for two-position 

processing based on separable upsampling according to one embodiment of the present invention. 
Figure 9 is a block diagram illustrating a system for generating a simulated high resolution image for two-position 
processing based on non-separable upsampling according to one embodiment of the present invention. 
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Figure 1 0 Is a block diagram illustrating a system for generating a simulated high resolution image for four-position 
processing according to one embodiment of the present invention. 

Figure 11 Is a block diagram illustrating the comparison of a simulated high resolution image and a desired high 
resolution image according to one embodiment of the present Invention. 

Figure 12 is a diagram illustrating the effect in the frequency domain of the upsampling of a sub-frame according 
to one embodiment of the present invention. 

Figure 13 is a diagram illustrating the effect in the frequency domain of the shifting of an upsampled sub-frame 
according to one embodiment of the present Invention. 

Figure 14 is a diagram illustrating regions of Influence for pixels In an upsampled image according to one embod- 
iment of the present Invention. 

Figure 15 is a diagram illustrating the generation of an initial simulated high resolution image based on an adaptive 
multi-pass algorithm according to one embodiment of the present invention. 

Figure 16 is a diagram illustrating the. generation of correction data based on an adaptive multi-pass algorithm 
according to one embodiment of the present invention. 

Figure 1 7 is a diagram illustrating the generation of updated sub-frames based on an adaptive multi-pass algorithm 
according to one embodiment of the present invention. 

Figure 18 is a diagram Illustrating the generation of correction data based on an adaptive multi-pass algorithm 
according to another embodiment of the present invention. 

20 Description of the Preferred Embodinfients 

[0007] In the following detailed description of the preferred embodiments, reference is made to the accompanying 
drawings, which form a part hereof, and in which is shown by way of illustration specific embodiments in which the 
invention may be practiced. It is to be understood that other embodiments may be utilized and structural or logical 
25 changes may be made without departing from the scope of the present invention. The following detailed description, 
therefore, is not to betaken in a limiting sense, and the scope of the present Invention is defined by the appended claims. 

I. Spatial and Temporal Shifting of Sub-frames 

30 [0008] Some display systems, such as some digital light projectors, may not have sufficient resolution to display 
some high resolution images. Such systems can be configured to give the appearance to the human eye of higher 
resolution images by displaying spatially and temporally shifted lower resolution Images. The lower resolution images 
are refen'ed to as sub-frames. A problem of sub-frame generation, which is addressed by embodiments of the present 
invention, is to determine appropriate values for the sub-frames so that the displayed sub-frames are close in appear- 

35 ance to how the high-resolution image from which the sub-frames were derived would appear if directly displayed. 
[0009] One embodiment of a display system that provides the appearance of enhanced resolution through temporal 
and spatial shifting of sub-frames is described In the above-cited U.S. patent applications, and is summarized below 
with reference to Figures 1-4E. 

[0010] Figure 1 is a block diagram illustrating an image display system 10 according to one embodiment of the 
40 present invention. Image display system 1 0 facilitates processing of an image 1 2 to create a displayed image 1 4. Image 
12 is defined to include any pictorial, graphical, and/or textural characters, symbols, illustrations, and/or other repre- 
sentation of infomriation. Image 12 Is represented, for example, by image data 16. Image data 16 includes individual 
picture elements or pixels of image 12. While one image is illustrated and described as being processed by image 
display system 10, it is understood that a plurality or series of images may be processed and displayed by image 
45 display system 10. 

[0011] In one embodiment, image display system 10 includes a frame rate conversion unit 20 and an Image frame 
buffer 22, an image processing unit 24, and a display device 26. As described below, frame rate conversion unit 20 
and image frame buffer 22 receive and buffer image data 1 6 for image 12 to create an image frame 28 for image 12. 
Image processing unit 24 processes image frame 28 to define one or more image sub-frames 30 for image frame 28, 

50 and display device 26 temporally and spatially displays image sub-frames 30 to produce displayed image 1 4. 

[0012] Image display system 1 0, including frame rate conversion unit 20 and/or image processing unit 24, includes 
hardware, software, firmware, or a combination of these. In one embodiment, one or more components of image display 
system 1 0, including frame rate conversion unit 20 and/or image processing unit 24, are included in a computer, com- 
puter server, or other microprocessor-based system capable of performing a sequence of logic operations. In addition, 

55 processing can be distributed throughout the system with individual portions being implemented in separate system 
components. 

[0013] Image data 16 may include digital image data 161 or analog image data 162. To process analog image data 
162, image display system 10 Includes an analog-to-digital (A/D) converter 32. As such, A/D converter 32 converts. 



3 



EP 1503 335 A1 



analog image data 162 to digital form for subsequent processing. Tlius, image display system 10 may receive and 
process digital image data 161 and/or analog image data 162 for image 12. 

[0014] Frame rate conversion unit 20 receives image data 16 for image 12 and buffers or stores image data 16 in 
image frame buffer 22. More specifically, frame rate conversion unit 20 receives Image data 1 6 representing individual 

5 lines or fields of image 12 and buffers image data 16 in image frame'buffer 22 to create image frame 28 for image 12. 
image frame buffer 22 buffers image data 16 by receiving and storing all of the image data for image frame 28, and 
frame rate conversion unit 20 creates Image frame 28 by subsequently retrieving or extracting all of the image data 
for image frame 28 from image frame buffer 22. As such, Image frame 28 is defined to include a plurality of individual 
lines or fields of image data 16 representing an entirety of image 12. Thus, image frame 28 includes a plurality of 

10 columns and a plurality of rows of individual pixels representing image 1 2. 

[0015] Frame rate conversion unit 20 and image frame buffer 22 can receive and process image data 16 as progres- 
sive Image data and/or Interlaced Image data. With progressive Image data, frame rate conversion unit 20 and image 
frame buffer 22 receive and store sequential fields of image data 1 6 for image 1 2. Thus, frame rate conversion unit 20 
creates image frame 28 by retrieving the sequential fields of image data 16 for image 12. With interlaced image data, 

15 frame rate conversion unit 20 and image frame buffer 22 receive and store odd fields and even fields of image data 
16 for image 12. For example, all of the odd fields of image data 16 are received and stored and all of the even fields 
of Image data 16 are received and stored. As such, frame rate conversion unit 20 de-Interlaces image data 16 and 
creates image frame 28 by retrieving the odd and even fields of image data 1 6 for image 12. 
[0016] Image frame buffer 22 includes memory for storing image data 16 for one or more image frames 28 of re- 

20 spectlve images 12. Thus, image frame buffer 22 constitutes a database of one or more image frames 28. Examples 
of image frame buffer 22 include non-volatile memory (e.g., a hard disk drive or other persistent storage device) and 
may include volatile memory (e.g., random access memory (RAM)). 

[0017] By receiving image data 16 at frame rate conversion unit 20 and buffering image data 16 with image frame 
buffer 22, input timing of image data 1 6 can be decoupled from a timing requirement of display device 26. More spe- 
25 cifically, since image data 1 6 for image franie 28 Is received and stored by image frame buffer 22, image data 1 6 can 
be received as Input at any rate. As such, the frame rate of Image frame 28 can be converted to the timing requirement 
of display device 26. Thus, image data 1 6 for Image frame 28 can be extracted from image frame buffer 22 at a frame 
rate of display device 26. 

[001 8] In one embodiment, image processing unit 24 includes a resolution adjustment unit 34 and a sub-frame gen- 
30 eratlon unit 36. As described below, resolution adjustment unit 34 receives Image data 16 for image frame 28 and 
adjusts a resolution of Image data 1 6 for display on display device 26, and sub-frame generation unit 36 generates a 
plurality of image sub-frames 30 for Image frame 28. More specifically, image processing unit 24 receives Image data 
16 for image frame 28 at an original resolution and processes Image data 16 to increase, decrease, and/or leave 
unaltered the resolution of image data 16. Accordingly, with image processing unit 24. image display system 10 can 
35 receive and display Image data 1 6 of varying resolutions. 

[0019] Sub-frame generation unit 36 receives and processes image data 1 6 for Image frame 28 to define a plurality 
of image sub-frames 30 for image frame 28. If resolution adjustment unit 34 has adjusted the resolution of Image data 
16, sub-frame generation unit 36 receives image data 16 at the adjusted resolution. The adjusted resolution of image 
data 16 may be Increased, decreased, or the same as the original resolution of image data 16 for image frame 28. 
40 Sub-frame generation unit 36 generates Image sub-frames 30 with a resolution which matches the resolution of display 
device 26. Image sub-frames 30 are each of an area equal to image frame 28. Sub-frames 30 each include a plurality 
of columns and a plurality of rows of individual pixels representing a subset of image data 16 of image 12, and have 
a resolution that matches the resolution of display device 26. 

[0020] Each image sub-frame 30 includes a matrix or array of pixels for image frame 28. Image sub-frames 30 are 
^ spatially offset from each other such that each Image sub-frame 30 includes different pixels and/or portions of pixels. 
As such, image sub-frames 30 are offset from each other by a vertical distance and/or a horizontal distance, as de- 
scribed below. 

[0021] Display device 26 receives image sub-frames 30 from image processing unit 24 and sequentially displays 
Image sub-frames 30 to create displayed Image 1 4. More specifically, as image sub-frames 30 are spatially offset from 
50 each other, display device 26 displays Image sub-frames 30 in different positions according to the spatial offset of 
image sub-frames 30, as described below. As such, display device 26 alternates between displaying image sub-frames 
30 for image frame 28 to create displayed image 14. Accordingly, display device 26 displays an entire sub-frame 30 
for image frame 28 at one time. 

[0022] In one embodiment, display device 26 perfomns one cycle of displaying image sub-frames 30 for each image 
55 frame 28. Display device 26 displays image sub-frames 30 so as to be spatially and temporally offset from each other. 
In one embodiment, display device 26 optically steers image sub-frames 30 to create displayed image 14. As such, 
individual pixels of display device 26 are addressed to multiple locations. 

[0023] In one embodiment, display device 26 Includes an image shifter 38. Image shifter 38 spatially alters or offsets 
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the position of image sub-frames 30 as displayed by display device 26. More specifically, image shifter 38 varies the 
position of display of image sub-frames 30, as described below, to produce displayed image 14. 
[0024] In one embodiment, display device 26 includes a light modulator for modulation of incident light. The light 
modulator includes, for example, a plurality of micro-mirror devices arranged to fonn an array of micro-mirror devices. 
5 As such, each micro-min-or device constitutes one cell or pixel of display device 26. Display device 26 may form part 
of a display, projector, or other imaging system. 

[0025] In one embodiment, image display system 10 includes a timing generator 40. Timing generator 40 commu- 
nicates, for example, with frame rate conversion unit 20, image processing unit 24, including resolution adjustment 
unit 34 and sub-frame generation unit 36, and display device 26, including Image shifter 38. As such, timing generator 
10 40 synchronizes buffering and conversion of image data 16 to create image frame 28, processing of image frame 28 
to adjust the resolution of image data 16 and generate Image sub-frames 30, and positioning and displaying of image 
sub-frames 30 to produce displayed image 1 4. Accordingly, timing generator 40 controls timing of image display system 
1 0 such that entire sub-frames of image 12 are temporally and spatially displayed by display device 26 as displayed 
image 1 4. 

15 [0026] In one embodiment, as illustrated in Figures 2A and 2B, image processing unit 24 defines two image sub- 
frames 30 for Image frame 28. More specifically, image processing unit 24 defines a first sub-frame 301 and a second 
sub-frame 302 for Image frame 28. As such, first sub-frame 301 and second sub-frame 302 each include a plurality of 
columns and a plurality of rows of individual pixels 18 of image data 16. Thus, first sub-frame 301 and second sub- 
frame 302 each constitute an image data an'ay or pixel matrix of a subset of image data 16. 

20 [0027] In one embodiment, as illustrated In Figure 2B, second sub-frame 302 is offset from first sub-frame 301 by a 
vertical distance 50 and a horizontal distance 52. As such, second sub-frame 302 is spatially offset from first sub-frame 
301 by a predetemriined distance. In one Illustrative embodiment, vertical distance 50 and horizontal distance 52 are 
each approximately one-half of one pixel. 

[0028] As illustrated in Figure 2C, display device 26 alternates between displaying first sub-frame 301 in a first position 
25 and displaying second sub-frame 302 in a second position spatially offset from the first position. More specifically, 
display device 26 shifts display of second sub-frame 302 relative to display of first sub-frame 301 by vertical distance 
50 and horizontal distance 52. As such, pixels of first sub-frame 301 overlap pixels of second sub-frame 302. In one 
embodiment, display device 26 perfomns one cycle of displaying first sub-frame 301 In the first position and displaying 
second sub-frame 302 in the second position for image frame 28. Thus, second sub-frame 302 is spatially and tem- 
30 porally displayed relative to first sub-frame 301 . The display of two temporally and spatially shifted sub-frames in this 
manner is referred to herein as two-position processing. 

[0029] In another embodiment, as illustrated in Figures 3A-3D, image processing unit 24 defines four image sub- 
frames 30 for image frame 28. More specifically, image processing unit 24 defines a first sub-frame 301 , a second sub- 
frame 302, a third sub-frame 303, and a fourth sub-frame 304 for image frame 28. As such, first sub-frame 301 , second 
35 sub-frame 302, third sub-frame 303, and fourth sub-frame 304 each Include a plurality of columns and a plurality of 
rows of individual pixels 1 8 of image data 1 6. 

[0030] In one embodiment, as Illustrated in Figures 3B-3D, second sub-frame 302 is offset from first sub-frame 301 
by a vertical distance 50 and a horizontal distance 52, third sub-frame 303 is offset from first sub-frame 301 by a 
horizontal distance 54, and fourth sub-frame 304 is offset from first sub-frame 301 by a vertical distance 56. As such, 

40 second sub-frame 302, third sub-frame 303, and fourth sub-frame 304 are each spatially offset from each other and 
spatially offset from first sub-frame 301 by a predetennined distance. In one illustrative embodiment, vertical distance 
50, horizontal distance 52, horizontal distance 54, and vertical distance 56 are each approximately one-half of one pixel. 
[0031] As illustrated schematically in Figure 3E, display device 26 alternates between displaying first sub-frame 301 
in a first position P^, displaying second sub-frame 302 in a second position P2 spatially offset from the first position, 

45 displaying third sub-frame 303 in a third position P3 spatially offset from the first position, and displaying fourth sub- 
frame 304 in a fourth position P4 spatially offset from the first position. More specifically, display device 26 shifts display 
of second sub-frame 302, third sub-frame 303, and fourth sub-frame 304 relative to first sub-frame 301 by the respective 
predetemriined distance. As such, pixels of first sub-frame 301 , second sub-frame 302, third sub-frame 303, and fourth 
sub-frame 304 overlap each other. 

50 [0032] In one embodiment, display device 26 perfonns one cycle of displaying first sub-frame 301 in the first position, 
displaying second sub-frame 302 in the second position, displaying third sub-frame 303 in the third position, and dis- 
playing fourth sub-frame 304 in the fourth position for image frame 28. Thus, second sub-frame 302, third sub-frame 
303, and fourth sub-frame 304 are spatially and temporally displayed relative to each other and relative to first sub- 
frame 301 . The display of four temporally and spatially shifted sub-frames in this manner Is refen-ed to herein as four- 

S5 position processing. 

[0033] Figures 4A-4E Illustrate one embodiment of completing one cycle of displaying a pixel 181 from first sub- 
frame 301 in the first position, displaying a pixel 182 from second sub-frame 302 In the second position, displaying a 
pixel 183 from third sub-frame 303 In the third posrtion, and displaying a pixel 184 from fourth sub-frame 304 in the 
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fourth position. IVIore specifically, Figure 4A Illustrates display of pixel 181 from first sub-frame 301 In the first position, 
Figure 4B Illustrates display of pixel 1 82 from second sub-frame 302 in the second position (with the first position being 
illustrated by dashed lines), Figure 4C illustrates display of pixel 1 83 from third sub-frame 303 in the third position (with 
the first position and the second position being Illustrated by dashed lines), Figure 4D illustrates display of pixel 184 

5 from fourth sub-frame 304 in the fourth position (with the first position, the second position, and the third position being 
Illustrated by dashed lines), and Figure 4E illustrates display of pixel 1 81 from first sub-frame 301 in the first position 
(with the second position, the third position, and the fourth position being illustrated by dashed lines). 
[0034] Sub-frame generation unit 36 (Figure 1) generates sub-frames 30 based on image data in image frame 28. 
It will be understood by a person of ordinary skill in the art that functions performed by sub-frame generation unit 36 

10 may be Implemented in hardware, software, firmware, or any combination thereof. The Implementation may be via a 
microprocessor, programmable logic device, or state machine. Components of the present invention may reside in 
software on one or more computer-readable mediums. The term computer-readable medium as used herein Is defined 
to include any kind of memory, volatile or non-volatile, such as floppy disks, hard disks, CD-ROMs, flash memory, read- 
only memory (ROM), and random access memory. 

15 [0035] In one form of the invention, sub-frames 30 have a lower resolution than image frame 28. Thus, sub-frames 
30 are also referred to herein as low resolution images 30, and image frame 28 Is also referred to herein as a high 
resolution image 28. It will be understood by persons of ordinary skill in the art that the temns low resolution and high 
resolution are used herein in a comparative fashion, and are not limited to any particular minimum or rhaximum number 
of pixels. In one embodiment, sub-frame generation unit 36 Is configured to generate sub-frames 30 based on one of 

20 five algorithms. These five algorithms are referred to herein as the following: (1) nearest neighbor; (2) bilinear; (3) 
spatial domain; (4) frequency domain; and (5) adaptive multi-pass. 

[0036] The nearest neighbor algorithm and the bilinear algorithm according to one form of the Invention generate 
sub-frames 30 by combining pixels from a high resolution image 28. The spatial domain algorithm and the frequency 
domain algorithm according to onefomn of the Invention generate sub-frames 30 based on the minimization of a global 

25 error metric that represents a difference between a simulated high resolution Image and a desired high resolution 
image 28. The adaptive multi-pass algorithm according to one fomn of the invention generates sub-frames 30 based 
on the minimization of a local error metric. In one embodiment, sub-frame generation unit 36 includes memory for 
storing a relationship between sub-frame values and high resolution image values, wherein the relationship is based 
on minimization of an error metric between the high resolution Image values and a simulated high resolution image 

30 that Is a function of the sub-frame values. Embodiments of each of these five algorithms are described below with 
reference to Figures 5-1 8. 

li. Nearest Neighbor 

35 [0037] Figure 5 is a diagram illustrating the generation of low resolution sub-frames 30A and 30B (collectively refenred 
to as sub-frames 30) from an original high resolution image 28 using a nearest neighbor algorithm according to one 
embodiment of the present invention. In the illustrated embodiment, high resolution image 28 includes four columns 
and four rows of pixels, for a total of sixteen pixels HI -HI 6. In one embodiment of the nearest neighbor algorithm, a 
first sub-frame 30A is generated by taking every other pixel in a first row of the high resolution image 28, skipping the 

^ second row of the high resolution image 28, taking every other pixel in the third row of the high resolution Image 28, 
and repeating this process throughout the high resolution Image 28. Thus, as shown in Figure 5, the first row of sub- 
frame 30A includes pixels HI and H3, and the second row of sub-frame 30A includes pixels H9 and H11 . In one fonii 
of the invention, a second sub-frame 30B is generated in the same manner as the first sub-frame 30A, but the process 
begins at a pixel H6 that Is shifted down one row and over one column from the first pixel H1 . Thus, as shown in Figure 

45 5, the first row of sub-frame 308 includes pixels H6 and H8, and the second row of sub-frame 308 includes pixels HI 4 
andH16. 

[0038] In one embodiment, the nearest neighbor algorithm Is implemented with a 2x2 filter with three filter coefficients 
of "0" and a fourth filter coefficient of " 1 " to generate a weighted sum of the pixel values from the high resolution image. 
Displaying sub-frames 30A and 308 using two-position processing as described above gives the appearance of a 
50 higher resolution image. The nearest neighbor algorithm is also applicable to four-position processing, and Is not limited 
to images having the number of pixels shown In Figure 5. 

III. Bilinear 

55 [0039] Figure 6 Is a diagram Illustrating the generation of low resolution sub-frames 30C and 30D (collectively referred 
to as sub-frames 30) from an original high resolution image 28 using a bilinear algorithm according to one embodiment 
of the present invention. In the illustrated embodiment, high resolution image 28 includes four columns and four rows 
of pixels, for a total of sixteen pixels HI -HI 6. Sub-frame 30C Includes two columns and two rows of pixels, for a total 
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of four pixels L1 -L4. And sub-frame 300 includes two columns and two rows of pixels, for a total of four pixels L5-L8. 
[0040] In one embodiment, the values for pixels L1-L8 in sub-frames 30C and SOD are generated from the pixel 
values H1-H16 of Image 2B based on the following Equations l-VIII; 

Equation I 
L1 =(4H1 +2H2+2H5)/8 

Equation II 
L2 := (4H3 + 2H4 +2H7) / 8 



Equation III 
L3=(4H9 + 2Hi0+2H13)/8 

Equation IV 
L4^(4H11 +2H12+2H15)/8 



Equation V 
L5=(4H6 + 2H2 +2H5)/8 



Equation VI 
L6=(4H8 + 2H4+2H7)/8 



Equation VII 
L7 = (4H14 + 2H10 +2H13) / 8 

^ Equation VIII 

L8 = (4H1 6 + 2H1 2 +2H1 5) / 8 

45 [0041] As can be seen from the above Equations l-VIII, the values of the pixels L1 -L4 in sub-frame 30C are influenced 
the most by the values of pixels H1, H3, H9, and H11, respectively, due to the multiplication by four. But the values for 
the pixels L1-L4 in sub-frame 30C are also influenced by the values of diagonal neighbors of pixels H1 , H3, H9, and 
H11 . Similarly, the values of the pixels L5-L8 in sub-frame 30D are influenced the most by the values of pixels H6, H8, 
H14, and H16, respectively, due to the multiplication by four. But the values forthe pixels L5-L8 in sub-frame 30D are 

50 also influenced by the values of diagonal neighbors of pixels H6, HB, H14, and HIS. 

[0042] In one embodiment, the bilinear algorithm is implemented with a 2x2 filter with one filter coefficient of "0" and 
three filter coefficients having a non-zero value (e.g., 4, 2, and 2) to generate a weighted sum of the pixel values from 
the high resolution image. In another embodiment, other values are used for the filter coefficients. Displaying sub- 
frames 30C and 30D using two-position processing as described above gives the appearance of a higher resolution 

55 image. The bilinear algorithm Is also applicable to four-position processing, and is not limited to images having the 
number of pixels shown in Figure 6. 

[0043] In one fomi of the nearest neighbor and bilinear algorithms, sub-frames 30 are generated based on a linear 
combination of pixel values from an original high resolution image as described above. In another embodiment, sub- 
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frames 30 are generated based on a non-linear combination of pixel values from an original high resolution image. For 
example, if the original high resolution image Is gamma-con-ected, appropriate non-linear combinations are used in 
one embodiment to undo the effect of the gamma curve. 

5 IV. Systems for Generating Simulated High Resolution Images 

[0044] Figures 7-10 Illustrate systems for generating simulated high resolution images. Based on these systems, 
spatial domain, frequency domain, and adaptive multi-pass aigorlthms for generating sub-frames are developed, as 
described in further detail below. 

10 [0045] Figure 7 is a block diagram illustrating a system 400 for generating a simulated high resolution image 412 
from two 4x4 pixel low resolution sub-frames 30E according to one embodiment of the present invention. System 400 
Includes upsampling stage 402, shifting stage 404, convolution stage 406, and summation stage 41 0. Sub-frames 30E 
are upsampled by upsampling stage 402 based on a sampling matrix, M, thereby generating upsampled images. The 
upsampled images are shifted by shifting stage 404 based on a spatial shifting matrix, S, thereby generating shifted 

IS upsampled Images. The shifted upsampled Images are convolved with an interpolating filter at convolution stage 406, 
thereby generating blocked images 408. In the illustrated embodiment, the interpolating filter is a 2x2 filter with filter 
coefficients of "1 and with the center of the convolution being the upper left position in the 2x2 matrix. The Interpolating 
filter simulates the superposition of low resolution sub-frames on a high resolution grid. The low resolution sub-frame 
pixel data is expanded so that the sub-frames can be represented on a high resolution grid. The interpolating filter fills 

20 In the missing pixel data produced by upsampling. The blocked images 408 are weighted and summed by summation 
block 410 to generate the 8x8 pixel simulated high resolution image 412. 

[0046] Figure 8 is a block diagram illustrating a system 500 for generating a simulated high resolution image 51 2 for 
two-position processing based on separable upsampling of two 4x4 pixel low resolution sub-frames 30F and 30G 
according to one embodiment of the present invention. System 500 Includes upsampling stages 502 and 514, shifting 

25 stage 518, convolution stages 506 and 522, summation stage 508, and multiplication stage 510. Sub-frame 30F is 
upsampled by a factor of two by upsampling stage 502, thereby generating an 8x8 pixel upsampled image 504. The 
dark pixels in upsampled image 504 represent the sixteen pixels from sub-frame 30F, and the light pixels in upsampled 
Image 504 represent zero values. Sub-frame 30G is upsampled by a factor of two by upsampling stage 514, thereby 
generating an 8x8 pixel upsampled image 516. The dark pixels in upsampled image 516 represent the sixteen pixels 

30 from sub-frame 30G, and the light pixels In upsampled Image 51 6 represent zero values. In one embodiment, upsam- 
pling stages 502 and 514 upsample sub-frames 30F and 30G, respectively, using a diagonal sampling matrix. 
[0047] The upsampled image 51 6 is shifted by shifting stage 51 8 based on a spatial shifting matrix, S, thereby gen- 
erating shifted upsampled image 520. In the illustrated embodiment, shifting stage 518 performs a one pixel diagonal 
shift. Images 504 and 520 are convolved with an interpolating filter at convolution stages 506 and 522, respectively, 

35 thereby generating blocked images. In the Illustrated embodiment, the interpolating filter at convolution stages 506 
and 522 is a 2x2 filter with filter coefficients of "1 and with the center of the convolution being the upper left position 
in the 2x2 matrix. The blocked images generated at convolution stages 506 and 522 are summed by summation block 
508, and multiplied by a factor of 0.5 at multiplication stage 510, to generate the 8x8 pixel simulated high resolution 
image 512. The image data is multiplied by a factor of 0.5 at multiplication stage 510 because, in one embodiment, 

40 each of the sub-frames 30F and 30G is displayed for only half of the time slot per period allotted to a color. In another 
embodiment, ratherthan multiplying by a factor of 0.5 at multiplbatlon stage 510, the filter coefficients of the Interpolating 
filter at stages 506 and 522 are reduced by a factor of 0.5. 

[0048] In one embodiment, as shown in Figure 8 and described above, the low resolution sub-frame data is repre- 
sented by two separate sub-frames 30F and 30G, whteh are separately upsampled based on a diagonal sampling 
45 matrix (I.e., separable upsampling). In another embodiment, as described below with reference to Figure 9, the low 
resolution sub-frame data is represented by a single sub-frame, which is upsampled based on a non-diagonal sampling 
matrix (i.e., non-separable upsampling). 

[0049] Figure 9 is a block diagram illustrating a system 600 for generating a simulated high resolution Image 61 0 for 
two-position processing based on non-separable upsampling of an 8x4 pixel low resolution sub-frame 30H according 

50 to one embodiment of the present Invention. System 600 includes quincunx upsampling stage 602, convolution stage 
606, and multiplication stage 608. Sub-frame 30H is upsampled by quincunx upsampling stage 602 based on a quincunx 
sampling matrix, Q, thereby generating upsampled image 604. The dari< pixels in upsampled image 604 represent the 
thirty-two pixels from sub-frame 30H, and the light pixels In upsampled Image 604 represent zero values. Sub-frame 
301-i includes pixel data for two 4x4 pixel sub-frames for two-position processing. The dark pixels In the first, third, fifth, 

55 and seventh rows of upsampled image 604 represent pixels for a first 4x4 pixel sub-frame, and the daric pixels in the 
second, fourth, sixth, and eighth rows of upsampled image 604 represent pixels for a second 4x4 pixel sub-frame. 
[0050] The upsampled image 604 is convolved with an interpolating filter at convolution stage 606, thereby generating 
a blocked image. In the Illustrated embodiment, the interpolating filter Is a 2x2 filter with filter coefficients of "1", and 
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with the center of the convolution being the upper left position In the 2x2 matrix. The blocked image generated by 
convolution stage 606 Is multiplied by a factor of 0.5 at multiplication stage 608, to generate the 8x8 pixel simulated 
high resolution image 610. 

[0051] Figure 1 0 Is a block diagram illustrating a system 700 for generating a simulated high resolution image 706 
5 for four-position processing based on sub-frame 301 according to one embodiment of the present invention. In the 
embodiment illustrated in Figure 10, sub-frame 301 is an 8x8 an-ay of pixels. Sub-frame 301 includes pixel data for four 
4x4 pixel sub-frames for four-position processing. Pixels A1 -A1 6 represent pixels for a first 4x4 pixel sub-frame, pixels 
B1-B16 represent pixels for a second 4x4 pixel sub-frame, pixels C1-C16 represent pixels for a third 4x4 pixel sub- 
frame, and pixels D1-D16 represent pixels for a fourth 4x4 pixel sub-frame. 
10 [0052] The sub-frame 301 is convolved with an interpolating filter at convolution stage 702, thereby generating a 
blocked image. In the illustrated embodiment, the interpolating filter is a 2x2 filter with filter coefficients of "1 and with 
the center of the convolution being the upper left position In the 2x2 matrix. The blocked image generated by convolution 
stage 702 is multiplied by a factor of 0.25 at multiplication stage 704, to generate the 8x8 pixel simulated high resolution 
image 706. The image data is multiplied by a factor of 0.25 at multiplication stage 704 because, in one embodiment, 
15 each of the four sub-frames represented by sub-frame 301 Is displayed for only one fourth of the time slot per period 
allotted to a color. In another embodiment, rather than multiplying by a factor of 0.25 at multiplication stage 704. the 
filter coefficients of the Interpolating filter are correspondingly reduced. 

V. Generation of Sub-frames Based on Error Minimization 



[0053] As described above, systems 400, 500, 600, and 700 generate simulated high resolution images 412, 512, 
610, and 706, respectively, based on low resolution sub-frames. If the sub-frames are optimal, the simulated high 
resolution image will be as close as possible to the original high resolution image 28. Various en-or metrics may be 
used to detemriine how close a simulated high resolution image is to an original high resolution image, including mean 

25 square error, weighted mean square error, as well as others. 

[0054] Figure 1 1 is a block diagram illustrating the comparison of a simulated high resolution image 41 2/51 2/61 01 706 
and a desired high resolution image 28 according to one embodiment of the present invention. A simulated high res- 
olution image 41 2, 51 2,610, or 706, is subtracted on a pixel-by-pixel basis from high resolution Image 28 at subtraction 
stage 802. In one embodiment, the resulting error image data is filtered by a human visual system (HVS) weighting 

30 filter (W) 804. 1 n one fomri of the Invention, HVS weighting filter 804 filters the en^or image data based on characteristics 
of the human visual system. In one embodiment, HVS weighting filter 804 reduces or eliminates high frequency errors. 
The mean squared error of the filtered data is then detemnined at stage 806 to provide a measure of how close the 
simulated high resolution image 412, 512, 610, or 706 is to the desired high resolution image 28. 
[0055] In one embodiment, systems 400, 500, 600, and 700 are represented mathematically in an error cost equation 

35 that measures the difference between a simulated high resolution image 412, 512, 610, or 706, and the original high 
resolution image 28. Optimal sub-frames are Identified by solving the en-or cost equation for the sub-frame data that 
provides the minimum error between the simulated high resolution image and the desired high resolution image. In 
one embodiment, globally optimum solutions are obtained in the spatial domain and in the frequency domain, and a 
locally optimum solution is obtained using an adaptive multi-pass algorithm. The spatial domain, frequency domain, 

40 and adaptive multi-pass algorithms are described in further detail below with reference to Figures 12-18. 

VI. Spatial Domain 

[0058] A spatial domain solution for generating optimal sub-frames according to one embodiment is described in the 
45 context of the system 600 shown in Figure 9. The system 600 shown in Figure 9 can be represented mathematically 
in an error cost function by the following Equation iX: 



20 



Equation IX 



50 




55 



where: 



^ optimal low resolution data for sub-frame 30H; 
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J = error cost function to be minimized; 

n and i<= indices for identifying liigh resolution pixel locations for Images 604 and 610; 
lQ(k) = image data from upsampled image 604 at location k; 
/(n - k) = filter coefficient of the interpolating filter at a position n - k; and 
s h(n) = Image data for desired high resolution Image 28 at location n. 

[0057] The summation of "lQ(k)f[n - k)" In Equation IX represents the convolution of the upsampled Image 604 and 
the Interpolating filter, f, performed at stage 606 in system 600. The filter operation Is performed by essentially sliding 
the lower right pixel of the 2x2 interpolating filter over each pixel of the upsampled image 604. The four pixels of the 
upsampled image 604 within the 2x2 Interpolating filter window are multiplied by the corresponding filter coefficient (i. 
e., "1 " in the illustrated embodiment). The results of the four multiplications are summed, and the value for the pixel of 
the upsampled image 604 corresponding to the lower right position of the interpolating fitter Is replaced by the sum of 
the four multiplication results. The high resolution data, h(n), from tine iiigh resolution image 28 is subtracted from the 
convolution value, lQ(k) f{n - k), to provide an error value. The summation of the squared error over all of the high 
resolution pixel locations provides a measure of the error to be minimized. 

[0058] An optimal spatial domain solution can be obtained by taking the derivative of Equation IX with respect to 
each of the low resolution pixels, and setting it equal to zero as shown in the following Equation X: 

^ Equation X 

4^ = 0, tee 
3/o(t) 

25 where: 

0 = the set of quincunx lattice points. 

[0059] Thus, as can be seen from Equation X, the derivative is taken only at the set of quincunx lattice points, which 
30 correspond to the dark pixels In upsampled Image 604 In Figure 9. Inserting the equation for J given In Equation IX 
Into Equation X, and taking the derivative as specified In Equation X, results In the following Equation XI: 

Equation XI 

35 

^Z^(k)C^(t-k) = 2;A(n)/(n-t), teG 

k n 

40 [0060] The symbol, Off, In Equation XI represents the auto>conrelatlon coefficients of the Interpolating filter, f, as 
defined by the following Equation XII: 

Equation XII 

45 
50 

[0061] Equation XI can be put into vector fomri as shown In the following Equation XIII: 

Equation XIII 
C/Q = h^. tG0 

where: 
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Cff= matrix of auto-corretation coefficients of the interpolating filter, f. 

vector representing the unknown image data for sub-frame 30H, as well as "don't care" data (i.e., the Image 
data corresponding to the light pixels in upsampled image 604); 

h/= vector representing a filtered version of the simulated high resolution image 610 using the interpolating filter, f. 

5 

[0062] Deleting the rows and columns corresponding to "don't care" data (i.e., the data that is not in the set of qun- 
incunx lattice points, &), results in the following Equation XIV: 



Equation XIV 

15 

where: 

\i = vector representing only the unknown Image data for sub - frame 30H. 

20 [0063] The above Equation XIV Is a sparse non-Toeplitz system representing a sparse system of linear equations. 
Since the matrix of auto-correlation coefficients is known, and the vector representing the filtered version'of the simu- 
lated high resolution Image 61 0 Is known. Equation XIV can be solved to detennine the optimal image data for sub- 
frame 30H. In one embodiment, sub-frame generation unit 36 is configured to solve Equation XIV to generate sub- 
frames 30. 

25 

VII. Frequency Domain 

[0064] A frequency domain solution for generating optimal sub-frames 30 according to one embodiment is described 
in the context of the system 500 shown in Figure 8. Before describing the frequency domain solution, a few properties 
30 of the fast fourier transfomi (FFT) that are applicable to the frequency domain solution are described with reference 
to Figures 12 and 13. 

[0065] Figure 12 is a diagram illustrating the effect In the frequency domain of the upsampling of a 4x4 pixel sub- 
frame 30J according to one embodiment of the present invention. As shown in Figure 12, sub-frame 30J is upsampled 
by a factor of two by upsampling stage 902 to generate an 8x8 pixel upsampled image 904. The dark pixels in upsampled 

35 image 904 represent the sixteen pixels from sub-frame 30J, and the light pixels in upsampled image 904 represent 
zero values. Taking the FFT of sub-frame 30J results in image (L) 906. Taking the FFT of upsampled image 904 results 
in image (LJj 908. Image (L^/) 908 Includes four 4x4 pixel portions, which are image portion (L^) 91 OA, image portion 
{L2) 91 OB, Image portion (L^ 91 OC, and image portion (^4) 91 OD. As shown in Figure 12, image portions 910A-910D 
are each the same as image 906 (i.e., = = L3 = L4 = L). 

40 [0066] Figure 1 3 is a diagram illustrating the effect in the frequency domain of the shifting of an 8x8 pixel upsampled 
sub-frame 904 according to one embodiment of the present invention. As shown in Figure 13, upsampled sub-frame 
904 is shifted by shifting stage 1002 to generate shifted image 1004. Taking the FFT of upsampled sub-frame 904 
results in image (LJ) 1006. Taking the FFT of shifted image 1004 results in image [L^S) 1008. Image (L^fS) 1008 
includes four 4x4 pixel portions, which are image portion (/.Sy) 1010A, image portion (LS^) 1010B, image portion {LS3I 

45 1010c, and image portion (LS4) 1010D. As shown in Figure 13, image 1008 is the same as image 1006 multiplied by 
a complex exponential, W[\.e., LijS= IV- Z.^), where denotes pointwise multiplication. The values for the complex 
exponential, W, are given by the following Equation XV: 

Equation XV 

50 

55 where: 

= row coordinate in the FFT domain; 
k2 = column coordinate in the FFT domain; 
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Ms number of columns In the image; and 
N = number of rows In the image. 

[0067] The system 500 shown In Figure 8 can be represented mathematically in an error cost function by the following 
5 Equation XVI: 

Equation XVI 



10 



(l;,l;)= argmin y = argniin 2[fXl. + W,^^^ 



15 



20 



25 



30 



where: 

(L*A» L*b) = vectors representing the optimal FFTs of sub-frames 30F and 30G, respectively, shown in Figure 8; 
J = en'oi' cost function to be minimized; 

/ = index Identifying FFT blocks that are averaged (e.g., for image 908 in Figure 12, four blocks are averaged, with 
/=1 corresponding to block 91 OA, i=2 corresponding to block 91 OB, i=3 corresponding to block 91 OC, and /=4 
corresponding to block 91 OD); 

F = matrix representing the FFT of the Interpolating filter, t 

Lf^ = vector representing the FFT of sub-frame 30F shown in Figure 8; 

Lg = vector representing the FFT of sub-frame 30G shown in Figure 8; 

W = matrix representing the FFT of the complex coefficient given by Equation XV; 

H = vector representing the FFT of the desired high resolution image 28. 

[0068] The superscript "H" in Equation XVI represents the Hennitian (i.e., X" is the Hermitian of X). The "hat" over 
the letters In Equation XVI indicates that those letters represent a diagonal matrix, as defined in the following Equation 
XVH: 



Equation XVII 



35 



40 



X=^diag{X]= 



(Xi 0 
0 X, 



0 
0 



0 
0 



0 X, 



0 
0 



0 0 X 



45 



[0069] Taking the derivative of Eq uation XVI with respect to the complex conjugate of La and setting It equal to zero 
results In the following Equation XVIII: 



Equation XVIII 



50 



O^A i i i 



55 [0070] Taking the derivative of Equation XVI with respect to the complex conjugate of Lq and setting it equal to zero 
results in the following Equation XIX: 
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Equation XIX 

^ - -2w,?.h, =0 

i A D 

[0071 ] The horizontal bar over the letters in Equations XVIII and XIX indicates that those letters represent a complex 
conjugate (i.e., A represents the complex conjugate of A). 

[0072] Solving Equations XVIII and XtX for and Lg results in the following Equations XX and XXI 



Equation XX 



Equation XXI 

L,=A-'(C^BL,) 



[0073] Equations XX and XXI may be implemented In the frequency domain using pseudo-inverse filtering. In one 
embodiment, sub-frame generation unit 36 is configured to generate sub-frames 30 based on Equations XX and XXI. 

VIII. Adaptive Multi-Pass 

[0074] An adaptive multi-pass algorithm for generating sub-frames 30 according to one embodiment uses past errors 
to update estimates for the sub-frame data, and provides fast convergence and low memory requirements. The adaptive 
multi-pass solution according to one embodiment is described in the context of the system 600 ishown in Figure 9. The 
system 600 shown in Figure 9 can be represented mathematically in an error cost function by the following Equation 
XXII: 



Equation XXII 



W(n) = |e('')(n)f =^|;/W(k)/(n-k)-/E(n)j 



where: 



n = index identifying the current iteration; 
J")(n) = error cost function at iteration n; 
e("){n) = square root of the error cost function, J("^{n)\ 

n and k = indices for identifying high resolution pixel locations in images 604 and 610; 

lQ(">(k) = image data from upsampled image 604 at location k; 

/(n - k) = filter coefficient of the interpolating filter at a position n - k; and 

h(n) = image data for desired high resolution image 28 at location n. 

[0075] As can be seen from Equation XXII, rather than minimizing a global spatial domain error by summing over 
the entire high resolution Image as shown in Equation IX above, a local spatial domain error, which is a function of n, 
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15 



35 



40 



is being minimized. 

[0076] A least mean squares (LMS) algorithm Is used in one embodiment to detemiine the update, which is repre- 
sented in the following Equation XXIII: 

Equation XXIII 



where: 



e = the set of quincunx lattice points (i.e., the dark pixels in upsampled Image 604 In Figure 9); and 

a = sharpening factor. 

[0077] Taking the derivative of Equation XXII provides the value for the derivative in Equation XXIII, which is given 
In the following Equation XXIV: 

20 Equation XXIV 

= 2^2^a '(k)/(n k) - A(n) j/(n - 1) 

25 



[0078] In one embodiment, a block-LMS algorithm using the average gradient over a "region of influence" is used 
to perfonn the update, as represented by the following Equation XXV: 

30 



Equation XXV 



a/W{n) 



where: 
Q ^ region of influence 

[0079] Figure 14 is a diagram illustrating regions of influence (Q) 1106 and 1108 for pixels in an upsampled image 
1 1 00 according to one embodiment of the present invention. Pixel 11 02 of image 1 1 00 con'esponds to a pixel for a first 
45 sub-frame, and pixel 11 04 of image 1 1 00 corresponds to a pixel for a second sub-frame. Region 1 1 06, which includes 
a 2x2 array of pixels with pixel 1102 in the upper left corner of the 2x2 array, is the region of influence for pixel 1102. 
Similarly, region 1108, which includes a 2x2 array of pixels with pixel 1104 in the upper left comer of the 2x2 array, is 
the region of influence for pixel 1104. 

[0080] Figure 15 is a diagram illustrating the generation of an initial simulated high resolution image 1208 based on 
50 an adaptive multi-pass algorithm according to one embodiment of the present Invention. An Initial set of low resolution 
sub-frames 30K-1 and 30L-1 are generated based on an original high resolution image 28. In the illustrated embodi- 
ment, the initial set of sub-frames 30K-1 and 30L-1 are generated using an embodiment of the nearest neighbor algo- 
rithm described above with reference to Figure 5, The sub-frames 30K-1 and 30L-1 are upsampled to generate up- 
sampled image 1202. The upsampled image 1202 is convolved with an interpolating filter 1204, thereby generating a 
55 blocked image, whrch is then multiplied by a factor of 0.5 to generate simulated high resolution image 1208. In the 
illustrated embodiment, the interpolating filter 1204 is a 2x2 filter with filter coefficients of "1", and with the center of 
the convolution being the upper left position in the 2x2 matrix. The lower right pixel 1206 of the interpolating filter 1204 
is positioned over each pixel in image 1202 to detennine the blocked value for that pixel position. As shown in Figure 
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15, the lower right pixel 1206 of the interpolating filter 1204 is positioned over the pixel in the third row and fourth 
column of image 1202, which has a value of "0". The blocked value for that pixel position is determined by multiplying 
the filter coefficients by the pixel values within the window of the filter 1204, and adding the results. Out-of-frame values 
are considered to be "0". For the illustrated embodiment, the blocked value for the pixel in the third row and fourth 
5 column of Image 1202 is given by the following Equation XXVt 



[0081] The value in Equation XXVI is then multiplied by the factor 0.5, and the result (I.e., 5) is the pixel value for 

the pixel 1210 in the third row and the fourth column of the initial simulated high resolution image 1208. 

[0082] After the initial simulated high resolution image 1208 is generated, correction data is generated. Figure 16 is 

15 a diagram illustrating the generation of correction data based on the adaptive multi-pass algorithm according to one 
embodiment of the present invention. As shown in Figure 16, the initial simulated high resolution image 1208 is sub- 
tracted from the original high resolution image 28 to generate an error Image 1302. Connection sub-frames 1312 and 
1314 are generated by averaging 2x2 blocks of pixels in error image 1302. For example, the pixel 1308 in the first 
column and first row of error image 1 302 has a region of influence 1 304. The pixel values within the region of influence 

20 1304 are averaged to generate a first correction value (i.e., 0.75). The first correction value is used for the pixel in the 
first column and the first row of correction sub-frame 1312. Similarly, the pixel 131 0 in the second column and second 
row of error image 1302 has a region of influence 1306. The pixel values within the region of Influence 1306 are 
averaged to generate a second correction value (I.e., 0.75). The second correction value is used for the pixel in the 
first column and the first row of correction sub-frame 1314. 

25 [0083] The correction value in the first row and second column of connection sub-frame 1 31 2 (i.e., 1 .38) is generated 
by essentially sliding the illustrated region of influence box 1304 two columns to the right and averaging those four 
pixels within the box 1 304. The correction value in the second row and first column of correction sub-frame 1 31 2 (I.e., 
0.50) is generated by essentially sliding the illustrated region of influence box 1304 two rows down and averaging 
those four pixels within the box 1304. The correction value in the second row and second column of correction sub- 

30 frame 1312 (i.e., 0.75) is generated by essentially sliding the illustrated region of influence box 1304 two columns to 
the right and two rows down and averaging those four pixels within the box 1304. 

[0084] The correction value in the first row and second column of con-ection sub-frame 1314 (i.e., 0.00) is generated 
by essentially sliding the Illustrated region of influence box 1306 two columns to the right and averaging those pixels 
within the box 1306. Out-of-frame values are considered to be "0". The correction value in the second row and first 

35 column of correction sub-frame 1314 (i.e., 0.38) is generated by essentially sliding the illustrated region of influence 
box 1306 two rows down and averaging those pixels within the box 1306. The correction value in the second row and 
second column of correction sub-frame 1314 (i.e., 0.00) is generated by essentially sliding the illustrated region of 
influence box 1306 two columns to the right and two rows down and averaging those four pixels within the box 1306. 
[0085] The correction sub-frames 1312 and 1314 are used to generate updated sub-frames. Figure 17 is a diagram 

40 illustrating the generation of updated sub-frames 30K-2 and 30L-2 based on the adaptive multi-pass algorithm accord- 
ing to one embodiment of the present invention. As shown in Figure 1 7, the updated sub-frame 30K-2 is generated by 
multiplying the correction sub-frame 1312 by the sharpening factor, a, and adding the initial sub-frame 30K-1. The 
updated sub-frame 30L-2 is generated by.multiplying the connection sub-frame 1314 by the sharpening factor, a, and 
adding the initial sub^frame 30L-1 . In the illustrated embodiment, the sharpening factor, a, is equal to 0.8. . 

45 [0086] In one embodiment, updated sub-frames 30K-2 and 30L-2 are used in the next iteration of the adaptive multi- 
pass algorithm to generate further updated sub-frames. Any desired number of iterations may be performed. After a 
number of iterations, the values for the sub-frames generated using the adaptive multi-pass algorithm converge to 
optimal values. In one embodiment, sub-frame generation unit 36 is configured to generate sub-frames 30 based on 
the adaptive multi-pass algorithm. 

50 [0087] The embodiment of the adaptive multi-pass algorithm described above with reference to Figures 1 5-1 7 is for 
two-position processing. For four-position processing, Equation XXIV becomes the following Equation XXVII: 



Equation XXVI 



10 



(1 xO) + (1 x5)-i-(1 x5) + (1 xO)= 10 



Equation XXVn 
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where: 



A'T) = low resolution data for the four sub-frames 30; 



5 [0088] And Equation XXlll becomes the following Equation XXVIII: 

Equation XXVIII 

10 /-i)(t)./")(t) + a^^ 



[0089] For four-position processing, there are four sub-frames, so the amount of low resolution data is the same as 
the amount of high resolution data. Each high resolution grid point contributes one error, and there is no need to average 

15 gradient update as represented in Equation XXV above. Rather, the error at a given location directly gives the update. 
[0090] As described above, in one embodiment, the adaptive multi-pass algorithm uses a least mean squares (LMS) 
technique to generate correction data. In another embodiment, the adaptive multi-pass algorithm uses a projection on 
a convex set (POCS) technique to generate correction data. The adaptive multi-pass solution based on the POCS 
technique according to one embodiment is described in the context of the system 600 shown In Figure 9. The system 

20 600 shown in Figure 9 can be represented mathematically in an error cost function by the following Equation XXIX: 



Equation XXIX 



|e(n)I = 

30 where: 

e{n) = error cost function; 
n and lc= indices identifying high resolution pixel locations;^ 
!^) = image data from upsampled image 604 at location k; 
35 /(n - k) = filter coefficient of the Interpolating filter at a position n - k; and 
h(n) =: image data for desired high resolution Image 28 at location n. 



5]/e(k)/(n-k)-/z(n) 



[0091] A constrained set for the POCS technique Is defined by the following Equation XXX: 



40 



45 



Equation XXX 

C(n) = J/e(n): 



'X'c(k)/(n-k)-/i(n)jj<;7j 



where: 



50 C(n) = constrained set that includes all sub-frame data from upsampled image 604 that is bounded by parameter, 

Ti; and 

r\ = error magnitude bound constraint. 



[0092] The sub-frame pixel values for the cun^ent iteration are detemnined based on the following Equation XXXI: 



55 
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Equation XXXI 



r"(t)= 



(l-;i)4"'(t)+A^^||j^ e(n')>Ti (tee) 
/W(t) ein') = tj 



15 where: 

n = index identifying the current iteration; 

X = relaxation parameter; and 

m = norm of the coefficients of the interpolating filter. 

20 

[0093] The symbol, n^, in Equation XXXI represents the location in the region of influence, SI, where the error is a 
maximum, and Is defined by the following Equation XXXII: 

25 Equation XXXII 

n' = argmaxjn e Q : \e{nj^ 



30 

[0094] Figure 1 8 is a diagram illustrating the generation of correction data based on the adaptive multi-pass algorithm 
using a POCS technique according to one embodiment of the present invention. I n one embodiment, an initial simulated 
high resolution image 1208 is generated in the same manner as described above with reference to Figure 1 5, and the 
initial simulated high resolution image 1208 is subtracted from the original high resolution Image 28 to generate an 
35 error image 1302. The Equation XXXi above Is then used to generate updated sub-frames 30K-3 and 30L-3 from the 
data in en-or image 1 302. For the illustrated embodiment, it is assumed that relaxation parameter, X, in Equation XXXI 
Is equal to 0.5, and the error magnitude bound constraint, "n, is equal to 1 . 

[0095] With the POCS technique, rather than averaging the pixel values within the region of influence to determine 
a connection value as described above with reference to Figure 16, the maximum error, e(n*), within the region of 
40 influence is identified. An updated pixel value is then generated using the appropriate fomnula from Equation XXXI, 
which will depend on whether the maximum ennor, e(n*), within. the region of influence is greater than 1 . less than 1 , 
or equal to 1 (since ti = 1 for this example). 

[0096] For example, the pixel in the first column and first row of error image 1302 has a region of influence 1304. 
The maximum error within this region of influence 13041s 1 (i.e., e(n*) = 1). Refemngto Equation XXXI, for the case 

45 where e(n*) - 1 , the updated pixel value is equal to the previous value for this pixel. Refening to Figure 1 5, the previous 
value for the pixel in the first column and the first row of sub-frame 30K-1 was 2, so this pixel remains with a value of 
2 In updated sub-frame 30K-3. The pixel In the second column and second row of error image 1302 has a region of 
Influence 1306. The maximum error within this region of influence 1306 is 1 .5 (i.e., e(n*) = 1 .5). Referring to Equation 
XXXI, for the case where e(n*) > 1 , the updated pixel value Is equal to half the previous value for this pixel, plus half 
. 50 of the quantity (e(n*) - 1), which is equal to 1 .25. Referring to Figure 15, the previous value for the pixel in the first 
column andthefirst row of sub-frame 30L-1 was 2, so the updated value for this pixel Is 1.25 in updated sub-frame 30L-3. 
[0097] The region of influence boxes 1 302 and 1 304 are essentially moved around the en-or image 1 302 in the same 
manner as described above with reference to Figure 16 to generate the remaining updated values in updated sub- 
frames 30K-3 and 30L-3 based on Equation XXXI. 

55 [0098] Although specific embodiments have been illustrated and described herein for purposes of description of the 
prefen-ed embodiment, It will be appreciated by those of ordinary skill in the art that a wide variety of alternate and/or 
equivalent implementations may be substituted for the specific embodiments shown and described without departing 
from the scope of the present invention. Those with siclll in the mechanical, electro-mechanical, electrical, and computer 
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arts will readily appreciate that the present invention may be implemented In a very wide variety of embodiments. This 
application is intended to cover any adaptations or variations of the preferred embodiments discussed herein. There- 
fore, it is manifestly intended that this invention be limited only by the claims and the equivalents thereof. 

/ 

5 

Claims 

1 . A method of displaying an Image (14) with a display device (26), the method comprising: 

10 receiving image data (1 6) for the image; 

generating a first sub-frame (30) and a second sub-frame (30) corresponding to the image data, the first and 
the second sub-frames generated based on minimization of an error between the Image data and a simulated 
image (412, 512, 610, or 706); and 

alternating between displaying the first sub-frame in a first position and displaying the second sub-franrie in a 
15 second position spatially offset from the first position. 

2. The method of claim 1 , and further comprising: 

generating a third sub-frame (30) and a fourth sub-frame (30) corresponding to the image data, the third and 
20 the fourth sub-frames generated based on minimization of the error between the Image data and the simulated 

Image; and 

wherein alternating between displaying the first sub-frame and displaying the second sub-frame further in- 
cludes alternating between displaying the first sub-frame in the first position, displaying the second sub-frame In 
25 the second position, displaying the third sub-frame in a third position spatially offset from the first position and the 

second position, and displaying the fourth sub-frame in a fourth position spatially offset from the first position, the 
second position, and the third position. 

3. The method of claim 2, wherein the simulated image is based on convolution of the first, the second, the third, and 
30 the fourth sub-frames with an interpolating filter (406, 506, 522, 606, or 702). 

4. The method of claim 3, wherein the Interpolating filter is a 2x2 interpolating filter with four filter coefficients, and 
wherein each of the filter coefficients has a value of one. 

35 5. The method pf claim 1 , wherein the simulated image is based on upsanripling of the first and the second sub- 
frames, thereby generating upsampled sub-frame data (504, 516, or 604). 

6. The method of claim 5, wherein the upsampled sub-frame data Includes first and second upsampled sub-frames 
(504 and 516), and wherein the simulated image is based on shifting of pixels in the first upsampled sub-frame, 

40 thereby generating a first shifted sub-frame (520), and wherein the simulated image is based on convolutions of 
the first shifted sub-frame and the second upsampled sub-frame with an interpolating filter (506 or 522), and where- 
in the simulated Image Is based on a summation of results of the convolutions. 

7. The method of claim 1 , wherein the error between the image data and the simulated innage represents a global 
^ error that is minimized in the spatial domain. 

8. The method of claim 1 , wherein the error between the image data and the simulated image represents a global 
error that is minimized in the frequency domain. 

50 9. The method of claim 1 , wherein the en^or between the image data and the simulated Image represents a local ennor 
that Is adaptively minimized In the spatial domain via a plurality of Iterations. 

10. A system (10) for displaying an image (14), the system comprising: 

55 a buffer (22) adapted to receive image data (1 6) for the image; 

an Image processing unit (24) configured to define first and second sub-frames (30) corresponding to the 
image data based on minimization of a difference between the image data and a simulated image (412, 512, 
610, or 706); and 
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a display device (26) adapted to alternately display the first sub-frame in a first position and the second sub- 
frame In a second position spatially offset from the first position. 
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